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Abstract. In the present work, we study correlations between stellar fundamental parameters and the oscillation 
amplitude for S Scuti stars . We present this study for a sample of 17 selected S Scuti stars belonging to 5 young 
open clusters. Taking advantage of properties of S Scuti stars in clusters, we correct the photometric parameters 
of our objects for the effect of fast rotation. We confirm the benefit of applying such corrections in this kind of 
studies. In addition, the technique used for this correction allows us to obtain an estimation of stellar parameters 
like the angle of inclination and the rotation rate, usually not accessible. A significant correlation between the 
parameter i! (estimation of the angle of inclination of the star) and the oscillation amplitude is found. A discussion 
and interpretation of these a priori surprising results is proposed, in terms of a possible selection rule for oscillation 
modes of S Scuti stars . 



Key words. Stars: 5 Set 
- Galaxies: star clusters 



Stars: rotation - Stars: statistics - Stars: oscillations - Stars: fundamental parameters 



1. Introduction 

The 5 Scuti stars are variables located in the lower part 
of the Cepheid instability strip with spectral types from 
A2 to FO. These pulsators, when belonging to the Main 
Sequence, seem particularly suitable for determining the 
extent of the convective core and internal rotation rate, 
and thereby for probing poorly understood hydrodynam- 
ical processes occurring in deep stellar interiors. Within 
the last decade great efforts have been made in develop 
ing the seismology of S Scuti stars (Breger 2000; Handler 
200C). However, several aspects of the pulsating behav- 
ior of these stars (e.g. the number and the determination 
of the excited modes) within the instability strip are not 
completely understood. 

Several statistical works investigated potential correla- 
tions between the oscillation amplitude and different stel- 
lar parameters. Particularly, a correlation between the am- 
plitude, period and absolute magnitude of low amplitude 
8 Scuti stars has been established (Antoncllo et al. 1981). 



Table 1. Main characteristics of the 5 open clusters used 
in this work. Different columns represent the cluster name, 
distance modulus, metallicity (columns 3 and 4), redden- 
ing and age. 



Cluster 


(to - M) v 


[M/H] 


Z 


E{B 2 - Vi) 


Age 




(mag) 






(mag) 


(Myr) 


Praesepe 


6.28 


0.170 


0.025 


0.0 


650 


Pleiades 


5.50 


-0.112 


0.014 


0.05 


130 


Hyades 


3.33 


0.143 


0.012 


0.0 


600 


Coma Ber 


4.70 


-0.048 


0.016 


0.0 


430 


a-Persei 


6.23 


-0.05 


0.016 


0.09 


90 



In a preliminary work for S Scuti stars in the Praesepe 



cluster, Suarez et al. (2001) have suggested the impor- 



tance of taking into account the effect of fast rotation as 
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proposed by Michel et al. (1999) in such studies. 
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Table 2. Indivual parallaxes for selected 5 Scuti stars of indexes are taken from the Rufener's cataloge (Rufcnci 



Hyades. it represents the parallax in miliarcseconds and 
(m — M) v is the resulting distance modulus. 



1999) 



Star 


7T 


(m 


-M) v 




(mas) 


(mag) 


HD 27397 


22.38 ± 0.37 


3.25 


±0.05 


HD 27459 


20.87 ±0.32 


3.40 


±0.10 


HD 27628 


22.30 ± 0.36 


3.25 


±0.05 



For a-Persei, the estimation of the global metallic- 
ity \M/H] by spectroscopi c measurements is taken from 
Pinsonneault et al. (1998| ). For the other four clusters, 



Here we extend this work to 5 open clusters, in or- 
der to estimate the impact of the mentioned technique in 
the studies of correlation. In addition, these tools allow 
us to obtain a rough estimation of several stellar parame- 
ters like the angle of inclination of the star or the rotation 
rate, which are not usually accessible. We present a study 
of single correlation between the parameters mentioned 
herebefore and the oscillation amplitude. For values indi- 
cating ia reasonable degree of corrclarion, we also search 



we used metallicities derived by M. Grenon from Geneva 
photometry of single stars with spectral type in the range 
F4-K3 ( |Grenon 2000[ |Robichon et al. 1999[ ). To correct the 
visual absorption due to interstellar matter, the relation 

A V = 3.6E(B 2 -Vi) 

was used. The intrinsic colours in the Geneva photometric 



sgree oi correlation, we 
for a possible functional dependence. 

This paper is organized as follows: in Section 2 we 
present the observational material and the main charac- 
teristics of the selected open clusters. Details on the com- 
putation of models are given in Section 3. In Section 4, 
we introduce the stellar parameters corrected for the ef- 
fect of Itast rotation. Section t> is dedicated to describe the 
technique ot correlation and the analysis ol the results. 
Finally, we bring up our conclusions in section 6. 



system were obtained from Nicolct (1981). 

The distance moduli for Coma (4.70 ± 0.04 mag) and 
Praesepe (6.28 ±0.12 mag) are derived from mean cluster 
parallaxes computed u sing Hipparcos intermediate data 
( Robichon et al. 1999 ). In the case of Pleiades, the esti- 
mation of the distance modulus is still a subject of dis- 
cussion. The Hipparcos values (5.36 ± 0.06) and the one 
obtained by Pinsonneault et al. (1998] ), 5.60 ± 0.0 6 mag 
illustrate the ran ge of values found in the literature ( Stella 
& Nissen 2001). For the Hyades, we have also access 



from literature, to the individual parallaxes of our tar- 
get 5 Scuti stars (dc Bruijne et al. 2001) which are 
listed in Table @. For the other stars we use an estima- 



tion of 3.33 ± 0.01 mag by Perryman et al. (199§| ) from 
Hipparcos observations. Finally for a-Persei, we used a 
distance modulus of 6.23 ± 0.06 from Pinsonneault et aT] 



(1998) 



The project ed velocity (Vs'mi) value s listed in Table || 
are taken from Uesugi fc Fukuda (1982). 



2. Observational data 



3. Stellar models and isochrones 



The oscillation amplitude (hereafter A osc ), as classically, is 
taken to be the highest peak-to-peak measurement in the 
observed light curve. Values for A osc used here are taken 



from Li fc Michel (1996 ) and references therein. The inter- 
pretation of this parameter is not straightforward; it can 
be considered as a maximum constructive interference of 
the observed modes. Thus, in such an interpretation, and 
for a given angle of inclination, it can be taken as an indi- 
cator of the energy of these modes. With such an estima- 
tion, it is evident that short time series might induce an 
underestimation of the amplitude, especially in presence 
of beating phenomena due to close frequencies. In order to 
take this effect into account, we associate with each am- 
plitude estimation a weight corresponding to the length 
of the time series. In Table | are listed such values where 
number 1 corresponds to effective time series shorter than 
4 hours; 2, between 4 and 20 hours; 3, between 20 and 
200 hours, and finally 4 corresponds to time series longer 
than 200 hours. The five young open clusters considered 
are: Praesepe, Pleiades, Hyades, Coma Ber and a-Persei, 
with ages between 90 and 700 Myr. Detailed information 
concerning metallicity, age, distance modulus and redden- 
ing are listed in Table |l|. Apparent magnitudes and colour 



The stellar evolutive models have been computed using 
the evolution code CESAM flMorel 1997|) . We consider 
input physics appropriate to the mass range covered by 
8 Scuti stars . Metallicity values (Z) used to compute 
models, are deduced from [M/H] assuming Y pl = 0.235 
and Z pr = as helium and heavy elements primordial 
concentrations, and an enrichment ratio of AY/ AZ = 2. 



We work with OPAL radiative opacity tables (Iglesias 
|fc Rogers 1996 ). For the atmosphere reconstruction, the 
Eddington's T(r) law is considered. Convection is de- 
scribed according to the classical mixing-length theory, 
with a solar calibrated mixing length parameter l m = 
1.61Hp, where Hp is the local pressure scale- height. For 
all the models, we take into account an overshooting of 
the mixed convective core over a distance a ov — 0.2Hp 



following the prescription of Schaller et al. (1992 ). 

Sets of evolutionary sequences were computed for 
masses between IMq and 5Mq from the zero-age Main 
Sequence to the subgiant branch. With representative se- 
quences for each cluster, we obtain the corresponding 
isochrones using the Geneva isochrone program. These 
isochrones were then transformed into a M v — (B2 — V\) 
diagram using the calibration of Kunzli et al. (1997) for 
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(B2 — V\ ) and the calibration of Schmidt-Kaler (1982| ) for 
M v . 

In Fig. [l] we present colour-magnitude (hereafter 
CM) diagrams with the observed data for each clus- 
ter. Continuous lines correspond to their associated 
isochrones, given the distances and ages listed in Table |l|. 
Potential binarity and the effect of fast rotation discussed 
in Section 4 are expected to induce systematic shifts to- 
ward higher luminosity and lower effective temperature, 
compared to single non-rotating stars. Thus, in order to 
compensate both effects, the fit of isochrones has been 
made adjusting them to the bottom envelope of the clus- 
ter in the CM diagrams. 

In the case of Pleiades, we use our best fit correspond- 
ing to a distance modulus of 5.5 mag, which is approxi- 
mately half way between the main sequence fitting result 
and Hipparcos distance (see Section 2). The corresponding 
age found is 130 Myr. 



Table 3. Correlation between parameters corrected for 
rotation and the A osc for our sample of 17 S Scuti stars . 
The Spearman Rank order coefficient is represented by r s ; 
r is the linear Pearson coefficient; a(r s ) er(r) are the errors 
in r s and r respectively, and p is a probability of deviation 
from null- hypothesis. 



Parameter 


r s 


a(r s ) 


p(r s ) 


r 


a(r) 


My 


0.4319 


0.1384 


0.0834 


0.4569 


0.1431 


B 2 -Vi 


0.2725 


0.0691 


0.2898 


0.2484 


0.0571 


•1 

i 


0.7481 


0.0447 


0.0005 


0.6553 


0.0447 


Vsmi 


0.7663 


0.0364 


0.0003 


0.7122 


0.0063 


V 


-0.2016 


0.0360 


0.4376 


-0.1810 


0.0290 


M/Mq 


0.3895 


0.1251 


0.1221 


0.4253 


0.1367 


J 


-0.1260 


0.0140 


0.629 


-0.1223 


0.0133 



4. Stellar parameters corrected for the effect of 
fast rotation 

We will not enter here in the details of the method, already 



described in Perez Hernandez et al. (1999). Following this 
paper, we apply the method to selected 5 Scuti stars for 
each of our five clusters. In Fig. ^ we illustrate such cor- 
rections for the five clusters. For a given Vsim, segments 
represent the potential position of the star, in absence of 
rotation, on the CM diagram. The position of this non- 
rotating copartners varies with the angle of inclination, i, 
and the rotation rate, lo. This last dimensionless quantity 
is defined by 

Q ~n 8GM 

with \l n — 



27R n 



where f2 is the angular rotational velocity of the star; fi c is 
the angular rotational velocity that a star with the same 



5.1. The errors 

We search for correlations between, on one side, the A osc 
parameter, and on the other side, stellar parameters es- 
sentially obtained after correcting for the effect of fast 
rotation. Here we discuss the error on these parameters. 

For A OSCy we estimate 0.005 and 0.01 as minimum and 
maximum typical errors and we calculate error bars as- 
suming a linear relation with the weights listed in Table ^[ 

Considering the stellar parameters, we have to deal 
with different error sources: input photometric in- 
dexes, Tasini values, distances, models and calibrations 
(isochrones), and of course, the method to correct param- 
eters for the effect of rotation. 

The photometric input parameters, basically the my 
and B2 — V1, have typical observational errors around 
10 -3 magnitudes. This is between 10 to 100 times smaller 



than the rotation effect found for these parameters (Perez 



mass (M) and polar radius (-ftp) would have if the centrifu- 
gal force balanced the gravitational attraction at the equa- 
tor. The non-rotating counterparts are thus, computed as 
a function of the angle of inclination, from i — 90° (w m i n ) 
to i m ; n obtained for ijj max = 0.95. 

For each selected 5 Scuti star we compute two associ- 
ated segments corresponding to corrections assuming an 
estimation of the error in Vsmi of ±10 %. Finally, mod- 
els at the intersection of segments with the isochrone 
will be our non-rotating counterparts, giving a set of cor- 
rected photometric parameters, like absolute magnitudes 
(M v ) and (B2 — V\) colour indices. In addition to this, the 
method gives us the possibility of obtaining an estimation 
of other stellar parameters like the angle of inclination of 
the star, the radial velocity, and the rotation rate, given 
in Table |. 

5. Searching for correlations with the oscillation 
amplitude 



Hernandez et al. 199E; Suarez et al. 2001). The error in 



the estimation of Tasini (estimated to 10%) is probably a 
dominant factor here. It can induce errors up to 0.15 and 
0.02 mag in M v and B2 — V\ corrected values respectively. 
The error on the fit of the isochrones (including those on 
the distance determination and on the calibration between 
B2 — Vi and T e ff) can be of the same order than the one 
associated to Vsmi. However they are expected to be sys- 
tematic, at least for objects within the same cluster. 

The errors associated more specifically to the mod- 
els and the correction for fast rotation are difficult to 
estimate, because we lack elements to compare with. 
Although for M v and B 2 — Vi, they are expected to be 
systematic to a good extent, for the rest of estimated stel- 
lar parameter we cannot guarantee it. In particular for i, 
V and uj, we can only pretend to obtain rough estimations, 
which will be indicated with a ('). 

We thus, consider that the error on Tasini is represen- 
tative of the non-systematic errors here, being the most 
important source of errors when searching for correlations. 
In addition, other observational aspects like unknown bi- 
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narity, differences with the mctallicity of the cluster and 
possible not detected differential extinction inside the clus- 
ter, can have an influence. 

For corrected parameters in Table |(| the error bars 
correspond to values obtained using ±10% in the observed 
T/sim. In Fig. ||we present the different diagrams with the 
parameters corrected for rotation versus A osc , with the 
error bars for both. 

5.2. The correlation technique 

We computed the Spearman rank-order correlation 
coefficient, r s . This coefficient is preferred to the classic 
Pearson one (also computed as a complement) as it is 
more robust to outliers and does not presuppose a linear 
relation. For the error made in the correlation coefficient, 
we use the following expression from standard statistics, 



cr(r s ) 



n - 1 



11 ry 
In 



where n is the dimension of the sample. In addition to 
this, we computed, for the Spearman's coefficient, the pa- 
rameter p which gives the probability of deviation from 
the null- hypothesis for this coefficient. 

On the other hand, for the r s values indicating a rea- 
sonable degree of correlation, we also search for a possible 
functional dependence. In such cases, we have used an 
orthogonal distance regression routine ( Press et al. 1989 , 
see) which combines data sets with errors in both vari- 
ables. 

5.3. Results and discussion 

We have applied the statistical methods described here- 
above to all parameters corrected for rotation computed 
in Section 4. In Table || are listed the coefficients obtained 
for correlations between these parameters and the A osc . 

We now analyze the results obtained for the absolute 
magnitude and the colour index. These two photometric 
parameters have been generally used in the literature for 
statistical works on S Scuti stars . Particularly, in a multi- 
correlate analysis made by Antonello et al. (1981), a signif- 
icant individual correlation between ^4 osc and M v is found. 
No strong correlation is found here neither for M v nor for 
B2 — Vt , however we notice that the correlation coefficient 
is approximately two times better for corrected values of 
M v than for non-corrected ones (r s = 0.26). They are of 
the same order in the case of B-2 — V\ (r s ~ 0.25). We thus 



conclude, following (Suarez et al. 2001), that it is worth- 
while to take into account such a correction in this kind 
of study. We attribute the low value obtained for the co- 
efficient in the case of the magnitude to the fact that our 
sample is still too limited. 

A less classical parameter for such studies is the pro- 
jected rotational velocity, Vsini, for which an unexpected 
significant correlation is found (r s = 0.76). This parameter 



Table 4. Linear functional dependence between Vsini, 
i'and ^4 sc- The linear fit parameters are a (xlO 3 ) and b 
following y = a^4 osc + b. Their corresponding errors are 
cr(a) (xlO 3 ) and <j(b). x 2 represents the minimum Chi- 
Square of the linear fit and q is a scalar between and 
1 giving the probability that a right model would give a 
value equal or larger than the observed \ 2 . 



y 


a 


a (a) 


(b) 


a(b) 


x 2 


q 


Vsmi 
i' 


1.88 
9.12 


1.12 
5.5 


-189.64 
-108.3 


305.430 
156.69 


3.997 
4.13 


0.997 
0.997 



has been considered before in |Li fc Michel (1999| ), however 
no correlation was calculated. This striking result has two 
possible explanations: a physical effect (linked to V' or 
maybe rather id') and/or a geometric effect (associated 
with i'). Making use of the different stellar parameters we 
have access to, we tried to go further to understand the 
origin of this correlation. We can notice that there is no 
significant correlation neither for the V parameter nor for 
the u)' parameter. 

On the other hand, for i', another significant correla- 
tion is found (r s = 0.75) such that, the higher is the value 
of A osc the larger is i'. These results suggest that the cor- 
relation observed for Tasini has a pure geometric origin. At 
this point, we tried to figure out if a possible bias could 
be induced by the method used to correct parameters for 
rotation. We can only figure out a possible influence of our 
limited sample of stars. In addition, the correlation found 
between V^sim and A osc does not involve the correction 
method. 

In the framework of the classical description of the os- 
cillations, which attributes to each mode a given spherical 
harmonic, the preceding result would first confirm that 
the variation of flux observed in 5 Scuti light curves are 
not dominated by the radial modes, which visibility coeffi- 
cients are independe nt of i' . This is in agreem ent with the 
results obtained by Hernandez et al. (1998 ) who showed 
that they could not find more than two radial modes per 
star, in oscillation spectra, for a group of selected stars in 
the Praesepe cluster. 

The A osc value being considered as the maximum con- 
structive interference of the observed oscillation modes, 
the present results would suggest that, at least, some 
modes with visibility coefficients increasing with i'are fa- 
vored in the S Scuti stars oscillations compared with other 
modes. 



As can be seen in Pesnell (1985), for instance, modes 
corresponding to this description are the sectoral modes 
(to = H), in opposition to the zonal modes ((£ — l,m = 
0),(t = 2, to = 0)) which visibility coefficient is low for i — > 
0, and (£ — 2, to = ±1) which visibility coefficient is low 
for i — > and for i — > 90. In this context, our results would 
suggest that the S Scuti oscillations are dominated, at least 
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to some extent, by modes with m = ±£. If confirmed, this 
interpretation would constitute an important new element 
in the problem of mode identification of S Scuti stars . 

However, we want to stress that the classical spherical 
harmonics system might not be the best context to inter- 
pret this result. It is well known now, that <5 Scuti stars are 
fast rotators (this is confirmed here where, for our sample 
of 5 Scuti stars , large V values are found) and that the 
effect of rotation affects their structure to a point where 
the description in terms of spherical harmonics might be- 
come questionable. An alternative representation is being 



developed by Lignieres (Lignieres 2001 ) by means of a non- 
perturbative numerical method. Their results obtained for 
homogeneous ellipsoids show modes with geometric hori- 
zontal features very different from classical spherical har- 
monics. These results are preliminary and have to be ap- 
plied to models closer to stellar internal structure, but 
they might be indicative of relevant geometric effects. 

For projected velocity and i', we have also searched 
for a functional dependence. In Table f| their linear fit 
parameters are listed. We have to keep in mind that the 
coefficients of these linear fits might be sensitive to even- 
tual systematic errors, which were neglected when search- 
ing for correlations (see Section 5). While good confidence 
values (q) are obtained, we can appreciate from values of 
X 2 , that scatter is important as well as errors in fit coef- 
ficients. This situation can be improved by increasing the 
number of target stars. 

Finally, no correlation is found for the mass. 

6. Conclusions 

For a sample of 17 5 Scuti stars , we report here the results 
of a search for correlations between the oscillation ampli- 
tude and a set of estimated fundamental parameters. 

In this process we confirm the interest of applying cor- 
rectio n for fast rotation, as suggested by ( Suarez et al 
200l| ) for this kind of study. 

We found a significant correlation coefficient for the 
projected velocity (Tasini) and show that it is due to the 
strong correlation between the oscillation amplitude and 
an estimation of the angle of inclination of the star, i' . 
We observe that the oscillation amplitudes increase with 
i' . This striking result is discussed in terms of the visibil- 
ity of modes. In the context of spherical harmonics, this 
suggests that the oscillations of 5 Scuti stars would follow 
a selection rule in favor of sectoral modes (i.e. m = ±£ 
modes). 

An alternative interpretation might turn up 
with future developments of the work initiated by 



Lignieres (Lignieres 2001). 

A future increase of our sample of S Scuti stars in clus- 
ters would help to refine the present results. 
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Table 5. Photometric data for the 17 6 Scuti stars of the 5 open clusters considered. Different columns represent 
the star HD (SAO) number, the absolute magnitude, the B 2 — V\ colour index, the oscillation amplitude A osc , the 
projected velocity Vsim,the spectral type and the class of luminosity, the time observation weight (see Section 2) and 
the cluster they belong to. 



Star 


M v 


B 2 - Vi 


^OSC 


Vsini 


Sp 


Weight 


Cluster 




(mag) 


(mag) 


(mag) 


(kms ) 








HD 74050 


1.62 


0.001 


0.015 


145 


A7V 


2 


Praesepe 


HD 74028 


1.67 


0.005 


0.020 


165 


A7III 


1 


Praesepe 


HD 73345 


1.87 


-0.001 


0.015 


95 


F0V 


1 


Praesepe 


HD 73175 


1.96 


0.037 


0.020 


170 


A9V 


2 


Praesepe 


HD 73798 


2.12 


0.063 


0.016 


170 


F0V 


4 


Praesepe 
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Fig. 1. Colour - Magnitude diagrams for the selected 5 open clusters. Filled circles and crosses represent respectively 
the observed S Scuti and other stars belonging to each cluster. The correction for distance and reddening is taken into 
account as mentioned in the text (see Section 2). 
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Fig. 2. Colour - Magnitude diagrams for the selected 5 open clusters. Filled circles represent the observed 5 Scuti . 
The continuous line is the isochrone corresponding to each cluster. Segments (two per star) represent the possible 
non-rotating counterparts for each observed star (see details in Section 4). 
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Table 6. Non-rotating associate models for the S Scuti stars in the 5 open clusters considered. The different columns 
represent the range of values obtained from the correction for the effect of rotation (see section 4) for: M v the absolute 
magnitude; B 2 — V\ the colour index; M the mass (in Mq); R the radius (in Rq); T e g the effective temperature (log, 
in K); L the luminosity (in Lq); uj' the rotation rate; v m f the rotation frequency (/iHz); i'the angle of inclination 
(deg), and finally g the gravity (log, in cms 2 ). 
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